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Abstract. We report magnetization (M) measurements close to the ferromagnetic quantum 
phase transition of the d-metal alloy Nii-^Va; at a vanadium concentration of x c « 11.4%. In 
the diluted regime (x > x c ), the temperature (T) and magnetic field (H) dependencies of the 
magnetization are characterized by nonuniversal power laws and display H/T scaling in a wide 
temperature and field range. The exponents vary strongly with x and follow the predictions of 
a quantum Griffiths phase. We also discuss the deviations and limits of the quantum Griffiths 
phase as well as the phase boundaries due to bulk and cluster physics. 



1. Introduction 

Magnetic quantum phase transitions (QPT) have been studied in transition metal alloys and 
in heavy-fermion compounds tuned, e.g., by pressure or chemical substitution. They still offer 
challenges to theory and experiment (see Ref. [I] for a recent review). Quantum critical behavior 
is signified by singularities in thermodynamic and transport properties. Usually, specific power 
laws with characteristic exponents have been predicted at the quantum critical point (QCP) 
for "clean" homogeneous systems, while "disordered" inhomogeneous systems, driven, e.g., by 
chemical substitution, may show different behavior [2]. In the case of metallic (itinerant) Heisen- 
berg magnets, a strong-disorder renormalization group [3] predicts an exotic infinite-randomness 
QCP accompanied by quantum Griffiths singularities |4j. At such a QCP, thermodynamic ob- 
servables are expected to be singular not just at criticality but in a finite region around the QCP 
called the quantum Griffiths phase (GP). This region features power laws (e.g., in the magnetic 
susceptibility, x ~ I 7 * -1 , and the magnetization, M ~ H x ) characterized by a nonuniversal 
Griffiths exponent A which varies with distance to the QCP. Quantum Griffiths singularities 
have attracted a lot of attention. Many heavy fermion compounds display anomalous power- 
laws in specific heat C(T) and x{T) [S]; and quantum Griffiths behavior was suggested as an 
explanation [B] . Recently, a more systematic variation of the exponents could be found at the 
ferromagnetic QPT of CePdi-^Rhj, [7]. 

To avoid additional complications due to the Kondo effect and to study a larger energy scale 
we recently investigated the simple fee transition metal alloy Nii-^Vx [5] as an example of an 
itinerant ferromagnet (FM) in which the transition temperature (T c = 630K for pure Ni) can be 
tuned to zero by chemical substitution. As explained by Friedel [9] the "disorder" is introduced 



because the charge contrast of the replacing vanadium atoms creates large defects yielding an 
inhomogeneous magnetization density. In contrast, diluting Ni with isoelectronic Pd does not 
lead to a strongly disordered scenario: Nii-^Pd^ remains ferromagnetic up to x c = 0.975 where 
it rather shows the signatures of a clean quantum critical point [10] . We showed in Ref. that 
magnetization and susceptibility above the critical vanadium concentration x c ~ 11.4% where 
T c is suppressed to indeed follow simple power laws with nonuniversal exponents that confirm 
the quantum Griffiths scenario over a wide temperature and magnetic field regime. At very 
low temperatures, deviations from the quantum Griffiths scenario hint at a cluster glass phase. 
Here, we provide additional details not shown in Ref. [11] . We demonstrate that i^/T-scaling 
holds for a wide concentration regime and show the scaling plots. In addition, we reveal how 
the impact of disorder is manifest in the original M data close to x c , and we show the details 
of the determination of the phase boundaries in order to better distinguish bulk behavior from 
individual cluster physics in this inhomogeneous system. 

Figure 1. Temperature (T) - concentration 
(x) phase diagram of Nii_ x V x showing the 
ferromagnetic (FM), paramagnetic (PM), 
quantum Griffiths (GP) and cluster glass 
(CG) phases. Closed and open diamonds 
mark T c and T c * as determined by linear and 
modified Arrott plots, respectively (data from 
[8] included). Triangles denote to T max from 
maxima in susceptibility. T cross is seen as 
lower limit of GP (see [11]). The dashed line 
is a linear extrapolation, the solid line is a 
tail fit describing the FM boundary, while the 
straight dotted line marks the onset of the 
CG. Inset shows the strong concentration x- 
variation of the exponent 7 from x(T), a — 1 
from M(H) and 7 from the H/T- scaling plot 
in the GP. 



2. Results 

Magnetization and ac-susceptibility measurements were performed on polycrystalline Nii_ z V z 
samples with x = 9 — 15% as described in Ref. [11] , An orbital contribution of Xorb = 
6 x 10~ 5 emu/mol has been subtracted from all data shown (M m = M — XorbH). 

Figure 1 shows the temperature-concentration phase diagram. For x < 11%, the critical 
temperature T c was determined by the standard Arrott analysis. Plots of M 2 vs. H/M as in 
Fig. 2 show straight parallel isotherms which implies M 2 = Mq (T) + cH/M as is common for 
itinerant magnets (in Fig. 2, only low T data are shown). T c is then extracted via the mean-field 
T-dependencies of Mq{T) and susceptibility (— c/Mq (T)) [12]. The resulting T c (x) can be simply 
extrapolated linearly (dashed line) from the high T c = 630K of nickel down to at x ~ 11% |8j. 

For x > 11%, the straightforward AP analysis does not longer work because the data in Fig. 
2 are not described by straight lines. Introducing "exponents" as in a classical critical regime 
leads to a "modified" Arrott plot [13J implying the behavior M 1 ^ = Mq /0 (T) + c(H/M) 1 ^. A 
good description for x > 11% of the M(H > 0.5T, T) data in a wide regime (outside any critical 
regime) can be achieved with (3 = 0.5 and 7(3;) < 1 [12J, as indicated by the dotted fit line in 




Fig. 2. The resulting transition temperatures T c * of these modified Arrott plots remain finite 
up to x = 15%, while other extrapolations in Fig. 2 using smaller H/M values would lead to 
smaller T c . 
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Figure 2. Probing Arrott plots (M 2 — 
Ml ~ {H/M) l h with 7 = 1) for various 
V-concentrations x. Far clarity only low 
temperature T = 2K data (with modified y- 
value) are shown. For x > 11% a modified 
Arrott plot (with 7(2) < 1) is a good 
description of the data above H ~ 0.5T as 
indicated by the dashed line. 
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Figure 3. Determination of transition tem- 
perature T max . (a) Differential susceptibility 
dM m /dH vs temperature T in small mag- 
netic fields H (symbol) displaying maxima 
at T max (H) for different x. (b) T max (H) vs 
magnetic field H. The dashed line indicates 
the linear extrapolation to determine T max = 
(H -> 0). (c) ac-susceptibility x'ac vs - 
temperature T in H = 0T with H ac = 0.1G 
at different frequencies v for x = 12.25%. 



In addition to Arrott plots, we analyze field-dependent maxima at T max (H) in the differ- 
ential susceptibility x(T") = dM(T)/dH indicating spin ordering or freezing as shown in Fig. 
3(a). Fig. 3(b) shows the linear extrapolation of T max (H) taken at 0.55T to 0.1T to deter- 
mine T max = T max {H — > 0). In particular for x = 12.25% a frequency dependent maximum 
at T max = 0.19K was determined by Xac{T,v = 380Hz) in zero field with H ac = 0.1G which 
hints at the onset of a cluster glass [11]. T max increases by 0.018K per decade in frequency 
|llj as shown in Fig. 3(c). Although a detailed study of the evolution with dilution x of the 
cluster growth and dynamics is still outstanding, we can already note the qualitative effects of 
disorder on the ferromagnetic ordered state for x > 11%. As is obvious in Fig. 1, the high 
and low field extrapolation lead to different transition temperatures (T c , T c * > T max ) hinting at 
cluster freezing for x > 11%. The ^-dependence of T max in the accessible temperature region 
is better described by an exponential (dotted line) rather than a power law. Also, a "tail" fit 
to (ln(T/T ) ~ (x c - x)-^, see [Hj) rather than a power law serves as a good description of 
the onset of FM order for data between about 9% and 11% leading to x c ~ 11.6% (solid line). 
The discrepancies between the various methods and the spin-glass like features at the lowest 
temperatures suggest that the real QCP is masked at very low T by ordering of clusters. 



Nonetheless, at sufficiently high temperatures (in the region T max < T < T c (0%)) cluster 
ordering does not seem to play a role, and various quantities display power laws. Figs. 4(a) and 
(b) present the H and T dependencies of the magnetization as M/H for various x. Fig. 4(b) 
shows essentially the susceptibility x, since x = M/H = dM/dH for low fields (H < 0.5T) 
and high T (T > 20K, T > T c ). While for x < 11% the negative slope in the log-log 
plot 7 = —dln(xm) / dln{T) increases with falling T towards T c , for x > 11%, x(T) follows 
a simple power law for 20K < T < 300A'. The exponent decreases from j(x = 11.4%) = 1 to 
7 ( x = 15%) = 0.04. M/H(H) follows a power law M/H ~ H a ~ l for high H. For x < 11%, 
where M(H) nearly saturates, the exponent 1 — a is close to 1, and therefore very different than 
7. However, for x > 12%, the high-field exponent I — a matches the susceptibility exponent 7. 
The deviations from a power law at low fields in Fig. 4(a) are due to the finite T limitations. 
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Figure 4. (a) Magnetic field H and (b) tem- 
perature T dependence of the magnetization 
M for a wide x regime. M/H(H, T) follows a 
power law (solid line) with the same exponent 
(a — 1) in (a) as (7) in (b) for all x > 11.6%. 




1 10 100 1000 1 10 100 1000 10 4 

H/T (G/K) H/T (G/K) 



Figure 5. H/T scaling plot showing some 
M(H, T) data within 10K-300K and 100G- 
50kG for different x. The line represents 
a fit using Y(z) (see text). Inset shows fit 
parameters fi SC ai and A vs x. 



Since both the M(T) and M(H) show power laws with the same exponent, simple H/T 
scaling is expected for x > 12%. Fig. 5 shows the scaling plot using the form M/H = 
H~^ ! 'Y([i sca [H /ksT) where Y is the scaling function and ^ sca i is a scaling moment for several 
x. All M(H,T) data for T > 14-fT collapse, confirming H/T scaling. The scaling function 
Y is well approximated by the form Y(z) = A' /(I + z -2 ) 7 / 2 where A' = A/p? 1 is a constant. 
This phenomenological form arises from simply combining the two limiting power laws with the 
same exponent 7 in the H — T plane, (M/H)- 1 = WY' 1 ~ [(n scal H) 2 + (k B T) 2 ]^/ 2 . Close 
to x = 11.6%, the quality of the collapse is less satisfactory. The resulting exponent 7 (which 
matches that obtained by a fit of x(T) for all x between 11.4% and 15%) is shown in the inset of 
Fig. 1. The scaling moment \x SC ai and amplitude A are shown in the inset of Fig. 5, demonstrating 
the growth of the typical cluster size and number with x — > x c . 

The consistent power laws, and in particular, the H/T scaling of M(H, T) are in excel- 
lent agreement with the predictions for a quantum Griffiths phase with Griffiths exponent 
X = a = 1 — 7. A critical concentration of x c = 11.4% can be identified from the condi- 
tion 7(x c ) = 1 (neglecting logarithmic terms). Fitting to power law 1 — 7(2;) ~ (x — x c ) v ^ as 
predicted by theory [31 yields x c = 11.6% with vtp = 0.42 as shown in the upper inset in Fig. 1. 
This value is in close agreement with the "tail" fit of T max 



3. Conclusions 

On the one hand, our results confirm that Nii-^V^ follows the scenario of an infinite-randomness 
QCP with a quantum Griffiths phase, as expected in an itinerant Heisenberg magnet [31 H]. The 
QCP at x c ~ 11.6% has been estimated by extrapolations from outside the critical region, where 
the cluster ordering is less disturbing (through "f(x c ) — > 1 and T max (x c ) — > 0). On the other 
hand, we see clear signs of cluster ordering towards x c , in particular deviations from scaling 
at lower temperatures (such as the upturns in Fig. 4(a) as well as model dependent transition 
temperatures for x > 11%. As discussed in Ref. [llj . the magnetization M(H,T > T max ) for 
x > 12% can be well described by an additional "Curie term" due to frozen clusters which 
exceeds the term due to the fluctuating (Griffiths) clusters below T cross (see Fig. 1). Such a 
change in low-temperature behavior was predicted to occur in itinerant Heisenberg systems due 
to the RKKY interactions [15] . A Griffiths phase with nonuniversal power laws at higher T (but 
below T c (0%)) combined with a cluster glass (CG) (indicated by maxima in x(^)) a t very low T 
has also been observed in other diluted compounds (CePdi-^Rhr [7J, URu9^Re-rSi9[T6]) close 
to a ferromagnetic transition with much lower T c and can be understood as a generic feature of 
this disordered itinerant QPT |14j . 
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